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Abstract 

A previous study of exactly solvable rationally-extended radial oscillator 
potentials and corresponding Laguerre exceptional orthogonal polynomials 
carried out in second-order supersymmetric quantum mechanics is extended 
to fcth-order one. The polynomial appearing in the potential denominator 
and its degree are determined. The first-order differential relations allowing 
one to obtain the associated exceptional orthogonal polynomials from those 
arising in a (/c — l)th-order analysis are established. Some nontrivial iden- 
tities connecting products of Laguerre polynomials are derived from shape 
invar iance. 
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1 Introduction 



Since the introduction of the first famihes Xi of Laguerre- and Jacobi-type rath- 
degree exceptional orthogonal polynomials (EOP), with n = 1, 2, 3, . . . , in the 
context of Sturm-Liouville theory [H |2] and the realization of their usefulness in 
building new exactly solvable rational extensions of known quantum potentials [3], 
a lot of work has been devoted to generalizing these families and the associated 
exactly solvable potentials, as well as to providing several different (but equivalent) 
approaches to the problem. 

After noting [B] that the first extended potentials enlarged the class of known 
translationally shape invariant potentials in supersymmetric quantum mechanics 
(SUSYQM) [U El Ej, it appeared convenient to use a SUSYQM technique to con- 
struct some additional examples of EOP and potentials [3 [8], in agreement with 
previous works on algebraic deformations of shape invariant potentials P| [TU]. 

Two distinct families of Laguerre and Jacobi Xm EOP, now labelled as type I 
and type II, were constructed for arbitrary large m pT|[T^ and their properties were 
thoroughly studied [131 E] • agreement with the general definition given in Refs. 
1 and 2, such nth-degree polynomials, with n = m, m + 1, m + 2, have the 
remarkable property of forming orthogonal and complete sets with respect to some 
positive-definite measure. They include as special cases the Xi EOP first considered 
in Refs. 1, 2, 3, and 7, as well as the X2 EOP introduced in Ref. 8. 

In addition, the X^ EOP were shown to be obtainable through several ap- 
proaches, such as the Darboux-Crum transformation fTbj \T6\ [T7] . the Darboux- 
Backlund one [TH], and the prepotential method |19] . 

Some novel procedures for generating translationally shape invariant potentials 
without recourse to EOP were also devised and led to examples of potentials |20l 
I2TI I22] , some of which were presented as new although they had already been built 
in Refs. 3 and 7. 

Recently, the X^ EOP (and associated potentials) were generalized to multi- 
indexed families Xmi,m2,...,mk^ obtained through the use of multi-step Darboux alge- 
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braic transformations [23], the Crum-Adler mechanism [23], higher-order SUSYQM 
[25] or multi-step Darboux-Backlund transformations [26] . 

In Ref. 25, in particular, the procedure for constructing rationally-extended ra- 
dial oscillator potentials and associated Laguerre EOP in higher-order SUSYQM was 
considered with special emphasis on second order. The number of distinct potentials 
and EOP families corresponding to a given degree fi of the polynomial arising in 
the potential denominator (which is some definite function of mi, m2, . . . , rrik) was 
determined for some small fi values and conjectured for higher ones. 

The purpose of the present work is to extend some of the results of Ref. 25 to 
order > 2, to establish relations between successive EOP, and to comment on 
shape invariance and some of its consequences. 

2 Potentials and EOP in First-order SUSYQM 

To start with, let us briefly review the construction of rationally-extended radial 
oscillator potentials and corresponding Laguerre EOP in first-order SUSYQM. 
As usual, the radial oscillator potential (in units h = 2m = 1) 

= ^c.V + (2.1) 

where u and / denote the oscillator frequency and the angular momentum quantum 
number, respectively, is defined on the half-line < a; < cxd. The corresponding 
Schrodinger equation has an infinite number of bound-state wavefunctions, which, 
up to some normalization factor, can be written as 

CX x'+^e-^-^'L?+^^(|u;x2) oc Vi{^)L^f{^), ly = 0,1,2, ... , (2.2) 

with 

z = ^ux', a = 1 + 1 r/K^) = -2^('"+'^e-^^, (2.3) 

and lI^\z) some Laguerre polynomial [27j. The associated bound-state energies are 
given by Ej)^ = uj{2v + / + §) = uj{2v + a + l). 
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In first-order SUSYQM [5] , one considers a pair of SUSY partners 

= = — ^ + (x) - E, H^-^ = AA^ = —^ + V^-^ (x) - E, 

dx^ dx'' ^2 4) 

= -^ + W(x), A = ^ + W(x), V^^\x) = W\x)tW'(x) + E, 
dx dx 

which intertwine with the first-order differential operators A and A^ as AH^'^^ = 

H^^^A and A'^H^~^ = H^'^^A'^. Here W{x) is the superpotential, which can be 

expressed as W{x) = — (log0(x))' in terms of a (nodeless) seed solution of the 

initial Schrodinger equation 



= E^{x), (2.5) 



E is the factorization energy, assumed smaller than or equal to the ground-state 
energy of V^+\ and a prime denotes a derivative with respect to x. Except in 
Sec. 5, we shall only consider here the case where E < Eq'^\ in which occurrence 
(f){x) is nonnormalizable, and we shall assume that the same holds true for (f)~^{x). 
Then H^^'^ and H^~'' turn out to be isospectral. 

For Vi{x), there are two different types of seed solutions (p{x) with such properties, 
namely 

0L(^) = x\iz)Ll:\-z) oc x'+^e^-^^L^+^Vlc^^'), (2-6) 
'Pllnix) = xnz)Ll;:'^\z) oc x-'e-'^-^'L^-'-^\lu;x% (2.7) 

with 

x]{z) = ^i(2«+i)e^, xYiz) = 2-i(2"-i)e-f-, (2.8) 

and corresponding energies Ej^ = —co{a + 2m + 1), Ej^ = —u{a — 2m — 1), respec- 
tively. Note that for type II, a must be greater than m. 

To obtain some rationally-extended radial oscillator potentials Vi^ext{x) with a 
given Z, we have to start from a conventional radial oscillator potential Vii{x) with 
some different V . We then get 

V^"-^ = Vv, V^~\x) = Vz,ext(x) +C = Vi{x) + Vi,,,,{x) + C, 



V^,rat(a;) = -2u { + 2z 



■ (a) ( -q) ^ 2 



(a) I (a) 
()m \ Qm 



(2.9) 



where a dot denotes a derivative with respect to z and C is some constant. According 
to the choice made for the seed function 0(a;), we may distinguish the two cases 

(z) r = /-l, = 0j_i,^, gt\z) = Lt-'^\-z), C = -u- (2.10) 
(u) /' = / + !, = 0j^i,^, g^::\z)=L^^--'\z), C = u. (2.11) 

The resulting extended potential is nonsingular on the half-line for any m = 1, 
2, . . . (and a large enough for type II). 

The bound-state wavefunctions ^^^-"^(x) of V*^^-* can be obtained by acting with 
A on those of V"*^+\ tl}\^^{x) oc rii\z)l}i \z) (with a' = /' + |), and are given by 

^^-\x)^-^^^y'f{z), n = m + v, = 0, 1, 2, . . . , (2.12) 

(2) 

where the nth-degree polynomial yli'\z) satisfies the differential equation 

yi^\z) = im-n)y^\z). 

(2.13) 

The orthonormality and completeness of tljlr'{x), = 0, 1, 2, . . . , on the half-line 
imply that the polynomials yi^\z), n = m + u, u = 0, 1, 2, . . . , form an orthogonal 
and complete set with respect to the positive-definite measure z'^e' -^{gi;:Y'dz. Ac- 
cording to the choice made for 0(x), such polynomials belong to the LI or L2 family 
and are denoted by L\^^{z) or L^^^^{z). They are conventionally normalized in 
such a way that their highest- degree term is given by (— z)"/[(r2 — m)!m!] multiplied 
by (-1)™ or 1, respectively [TT|[T^[Tl]f 

3 Potentials and EOP in fcth-order SUSYQM 

Going from first- to fcth-order SUSYQM [28l |29l [301 |3ll ESI |33] is achieved by 
replacing the first-order differential operators A, A^ by some fcth-order ones A^ 
in such a way that the SUSY partner Hamiltonians 

H^'^ = -^,+ V^'\x), if(^) = -^,+ V^'\^) (3.1) 

^It is worth noting, however, that in Ref. 15, there is an additional lactor — (a + n — 2to + 1) 
in the definition of ni-^)- 
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^ + a + 1 




•{") •■(") 

X g-m gm 
a)^^ + z- 



(a) 

gm 



(a) 
9m 



intertwine with them as AH^^^ = H^'^^A and A^H'^'^^ = H^^^A^. Here, we restrict 
ourselves to the reducible case, which could be formulated in the fcth-order parasu- 
persymmetric language as well [34j. The operator A can then be factorized into a 
product A^''^ A^'^"^^ . . . A^'^^ of first-order differential operators A^*^ = t" + W^^*^ 



X) 



i = 1, 2, . . . , k (and similarly for A'^). 

In terms of k seed functions 01, 02, • • • , 0fc of the starting Hamiltonian H^^\ the 
k functions W^^\x), i = 1, 2, . . . , k, can be expressed as W^^^x) = — (log(/)*^*)(a;))', 
where 

0«(x) = . . . A(^)0. = (3.2) 

W(01,02,...,0i-l) 

and W(0i, 02, • • • , 4>i) denotes the Wronskian of 0i(x), 02(a^), • • • , The po- 

tentials of the two partner Hamiltonians are linked by the relationship 

V('\x) = V^'\x) - 2^ log W(0i, 02, . . . 0fc). (3.3) 

As in Sec. 2, we start from some conventional radial oscillator potential V^^\x) = 
Vii{x) to get as a partner an extended potential with a given /, V^'^^{x) = Vi^ext{x) + 
C = Vi{x) + Vi^ratix) + C, up to somc additive constant C. Since we have two 
types of seed functions at our disposal and their order is irrelevant as far as the 
final potential is concerned, there are altogether k + 1 possibilities for the set of k 
seed functions, which we may denote by PII^^'', where q runs over 0, 1, . . . , /c. On 
assuming that the first q functions 0^ are of type I, we may write 

'xl{zW^?{-z) if ^ = 1,2, 
x]}{z)Lkr'\z) if^ = q + l,q + 2,...,k, 

with a' = /' + |, < mi < m2 < ■ • • < rrig, and < niq+i < mq^2 < ■ ■ < TUk- For 
a' large enough, these functions are nodeless on the half-line. 

The two pure cases l'^ and II'^, corresponding to q = k and q = 0, respectively, 
are very easy to deal with. For the former, on assuming V = I — k and using some 
simple properties of Wronskians [52] , we indeed obtain 

>V(0i,02,...,0.) = {^xf'''~'^)\xUT9'~^\z). 

(3.5) 

g^^\z) = W{Ll:[\-z), L^^:\-z), L^^:\-z)) . 
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Similarly, for the latter and I' = I + k, we get 

(3.6) 

Here W denotes Wronskians with respect to the variable z, while fi is the degree of 
the polynomial g'f'\z): 

gl^\z) = Cj^^z^ + lower order terms. (3.7) 

The values of n and C^""* can be easily found by replacing all Laguerre polynomials 
by their highest- degree term, e.g. Lm\z) by (— 2;)"Ym!, in the two Wronskians and 
evaluating the two resulting determinants. In both cases, fi = J2^=i " ^^{k — 1) 
and Cj^^ = ( — l)°"A(mi, m2, . . . , m^) / (mi!m2! . . . rrikl), where A(mi, 1712, ■ ■ ■ , rrik) = 
11^=4+1 ("^i ~ ^i) is ^ Vandermonde determinant of order k. In case l'^, a = 0, 
while in case II'^, a = mi + m2 + • • ■ + m^. 

In the mixed cases PIl'^"'', < g < A;, let us assume I' = I + k — 2q. On taking 
successively into account that 



x\}{z)=x'Az)z-'^'e-^, (3.8) 



and [27] 



d 



^Ll:'\-z) = L 
az 



(3.9) 



dz 

the Wronskian can be written as [26 



X W(L(:;)(-z), . . . , <)(-^), z-"'e-^LL-;;)(^), . . . , ^-"'e-^L(„7')(^)) 
= (a;x)'=(^-i)/2 (xJO ' (^-"'-'^+^e-^) det T^f ) , (3.10) 

elements of the kxk matrix F^"'' 

mer's symbol by 



where the matrix elements of the kxk matrix F/T are defined in terms of Pochham- 



^ + W-^4T;.^^^^W ifj = g + l,g + 2,...,fc, ^'-''^ 



for any i = 1, 2, . . . , k. 

To put W(0i, (f)2, ■ ■ ■ , (pk) into a form that interpolates between those for the pure 
cases, given in (13 .Sp and (I3.6p . it is convenient to factorize it as 



9, 



[a) 



{z) = z-^k-g)ik-,-i) detr(,°), 



(3.12) 



by using Eq. (13 .Sp again. Determining /z and C^"\ as defined in (13. 7p . is, however, 
not straightforward anymore, because the replacement of the Laguerre polynomials 
in (13. lip by their highest- degree term leads to a matrix with k — q proportional 
columns, whose determinant vanishes. 

To avoid such a drawback, it is appropriate to transform T^/^^ into another matrix 
f with equal determinant, for which such a problem does not occur. For any j = 1, 
2, . . . , k, the matrix elements of the latter may be defined by 

I'fH^ _ J if i = 1, 2, . . . ,g + 1, (QTQ\ 

^ ' 'EU.if-!;^)(rr^. iU = g + 2,q + 3,...,k. ^'-''^ 

The summation on the right-hand side of this equation can be easily carried out by 
using definition (13. lip and the equations |2Z] 

zLl:+'\z) + (m + l)4ti(^) = (m + a + l)Lt\z) 



L^J^-'Hz) + L^::Uz) = L^^^ 



(3.14) 



z . 



It can indeed be shown by induction over i running over g + 2, g + 3, . . . , A; that 

{L%tl~^\-z) if j = 1,2, ...,g, 

{m, + l),{m, -a'-t + q + 2),_,_i (3.15) 
xz'^-^Li-,;i^'^'\z) if J = g + 1, g + 2, . . . , fc. 

For gl^\z) = z~^^~''-'>'^^~'^~^'> det f with Laguerre polynomials replaced by their 
highest- degree term, we now directly get 

^ 1 1 
^ = ^mi- -q{q - 1) - -{k - q){k - q - 1) + q{k - g), (3.16) 

i=l 

- / 1 "^2, ■ ■ ■ , mg)A(mg+i, mg+2, ■ • • , m^) 

~ \ ^) III ' 

milm2\ . . . mk'. 

k (3.17) 
cr = mi + q{k- g), 

i=q+l 



where Eq. f l3.16p agrees with a result stated without proof in Ref. 24. As it is clear 
that Eqs. f l3.16p and fl3.17p . as well as the expression for the Wronskian in fl3.12p . 
include the corresponding results for pure cases (provided A(mi,m2, . . . , m^) = 1 
for g = 0), we may treat all q values simultaneously. 

It is then straightforward to see from Eq. fl3.3p that the two partner potentials 
V'^^^x) and V^'^\x) assume a form similar to that of V^^\x) and V^~\x) in Eq. 
(12. 9p . where we only have to substitute /i for m, I + k — 2q for and {k — 2q)uj for 
C. Provided l/^^) 

(x) is nonsingular, its bound-state wavefunctions ip^\x) satisfy 
Eqs. (I2.12P and (12.130 with fi substituted for m. We will now denote the orthogonal 
and complete EOP yii'\z), n = /i + z/, i/ = 0, 1, 2, . . . , (with respect to the 
positive-definite measure z^e'^ [gj^i'^) '^dz) by La,mifn2,...,nn:,n{z) , where mi, m2, ■ ■ ■ , 
niq belong to type I and mg+i, • • • , rnk to type II. We choose to normalize 

them so that their highest- degree term is (— z)"/[(r2— yu)!mi!m2! . . . m^!]. This agrees 
with the choice made in Ref. 23 for k = 2, q = 0, but differs from the conventional 
one for L\^^{z) (see at the end of Sec. 2). 

4 Relations between Successive Laguerre EOP 

The purpose of the present Section is to derive the first-order differential relations 
expressing the Laguerre EOP obtained in fcth-order SUSYQM in terms of those 
arising after k — 1 steps with the use of the first k — 1 seed functions 0i, 02, . . . , 
(t>k-i- Since the additional seed function (pk is of type I in the pure case l'^, but of 
type II in all the remaining cases, we have to treat them separately. 

In the former case, where we start from VJ_jt(x), we arrive after k — 1 steps at 
an extended potential Vi^i^extix) , which assumes a form similar to that of Vi^cxti^) 
but with gl^\z) replaced by gl^ where /ii = Yl^=i frii — ^{k — l){k — 2) . From 

the bound-state wavefunctions rii^i{z)y^n[ ^\z)/gl^ f^i = f^i, yUi + 1, of 

V5-i,ext(a;), we can obtain those of VJ,ext(a;), r]i{z)y^\z) / g''^\z) , n = /x, /i 1, . . . , 
by acting with the differential operator A^'^^ defined in Sec. 3, namely 




(4.1) 
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On using fl2.3p . we get the following equations 



/il ,^ a—l,mi 



mi,m2,.-->"ifc— ii'^l 
1 



a, mi ,m2,...,mf.,n V^/ ' 



(4.2) 



(4.3) 



where rii = ^i + i', n = fi + u, and any u = 0, 1, 2, . . . . The constant D appearing in 
(14.21) can be derived from the normahzation coefficients of the two wavefunctions and 
SUSYQM theory. It can ahernatively be obtained by comparing the highest-degree 
terms on both sides of (14. 2p . The result reads 

k-l 

D = {-ir^->^l[{mk-m,). (4.4) 



i=l 



In the latter case, where we start from Vi^k-2q{x) , < q < k,we arrive after k — 1 
steps at an extended potential V;+i^cxt(a;), associated with a polynomial glj^'^^\z), 
where /ii = ^Jl/ nii - |g(g - 1) - |(A; - g - 1)(A; - g - 2) + g(A; - g - 1). As 
in the previous case, from its bound-state wavefunctions r]i+i{z)ylil~^^\z) / gjj^^^\z) , 
ni = yUi, /ii + 1, . . . , those of V/^ext(2^) are obtained by application of 



A(k) I d 1 2« + 1 gi^ 

A^"' =LOX \ — -\ ^ 

dz 2 Az ^(°) 



(4.5) 



The resulting equation is now 



(4.6) 



with 



D 



(a+l) 

9 k 



(a) 



z— + a + 1 
dz 



^9. 



k-l 



-l)''~'^~^(a + iy + k-2q- rrik) J]^ (m^. - m^) 

i=q+l 



(4.7) 
(4.8) 
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5 Shape invariance and some of its consequences 

It is easy to see that if a rationally-extended radial oscillator potential Vi^cxti^) 
of the type considered in Sees. 3 and 4 has a ground-state wavefunction equal to 
'r]i{z)gl^'^^\z)/gj^\z) up to some constant multiplicative factor, then it is (transla- 
tionally) shape invariant. Such a property is based on the fact that if we consider 
a pair of partner potentials {V^'^\x),V^~\x)) with V^'^\x) = V;^cxt(a^) and if we 
take the above-mentioned wavefunction as factorization function and its eigenvalue 
Eq^ as factorization energy, then the corresponding superpotential W{x) can be 
separated into W{x) = Wi{x) + W2{x) with 

W^{x) = -ux - W,{x) = -ux[%^-%^], (5.1) 

and the partner potential is such that 

V^-\x) - if = [V,+i,ext(x) - i^r^T + 2uj. (5.2) 

Conversely, if Eq. f l5.2p is satisfied, then the lowest-degree EOP yli'\z) is propor- 
tional to gl^'~^^\z). 

Since the shape invariance of V;^ext(a;) has been explicitly shown in Ref. 26, in all 
cases < g < /c, we may write 



, , , ('_i)ELi'">-^9(9-i)-|('=-9){fc-'/-i) , 

a,m„m,,...,m,A ) A(mi , , . . . , m,) A(m,+i , m,+2 , • • • , ^fc) ^ ^ • ^ 



where the multiplicative factor has been found by comparing the highest-degree 
terms on both sides of the equation. 

We can now combine this result with Eqs. (Oj) . dO]) . flOj) . (Ojl . IK7} . and 
(14. Sp . where we now assume z/ = 0, rii = /ii, n = /i, to derive the two identities 



(a) 



9ii^ = -9t'^9lr'''\ (5.4) 
9lT"^ = {a + k-2q- m,)g^^^^'\^^^'\ (5.5) 
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valid ioi q = k and < g < fc, respectively. These are nontrivial (/ii + /x)th-degree 
relations connecting products of Laguerre polynomials. They can be easily checked 
for k = 1 and k = 2. In the former case, for instance, with gjj^^ = 1 and gjj^^ = gm\ 
as given in f l2.10p or (12. lip , the identities reduce to —Lm — Lm = 

-Lt\-z) and (« + 1)lL"-')(z) - = (a + 1 - m)Li7°-')(z) for g = 1 

and g = 0, respectively. 

6 Conclusion 

In the present work, we have determined the isospectral rationally-extended radial 
oscillator potentials resulting from conventional ones in fcth-order SUSYQM. We 
have found, in particular, the polynomial g^^\z) appearing in their denominator, as 
well as its degree /x. We have also established the first-order differential relations 
allowing one to obtain the corresponding Laguerre EOP from those arising in {k — 
l)th-order SUSYQM. Finally, we have proved some nontrivial identities connecting 
products of Laguerre polynomials, which are direct consequences of the extended 
potential shape invariance. 

A similar study could be carried out for Jacobi EOP and associated rationally- 
extended potentials. We hope to come back to this problem in a forthcoming work. 
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